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We consider systems of interacting spins and study the entanglement that can be localized, on 
average, between two separated spins by performing local measurements on the remaining spins. 
This concept of Localizable Entanglement (LE) leads naturally to notions like entanglement length 
and entanglement fluctuations. For both spin-1/2 and spin-1 systems we prove that the LE of a 
pure quantum state can be lower bounded by connected correlation functions. We further propose 
a scheme, based on matrix-product states and the Monte Carlo method, to efRciently calculate the 
LE for quantum states of a large number of spins. The virtues of LE are illustrated for various 
spin models. In particular, characteristic features of a quantum phase transition such as a diverging 
entanglement length can be observed. We also give examples for pure quantum states exhibiting 
a diverging entanglement length but finite correlation length. We have numerical evidence that 
the ground state of the antiferromagnetic spin-1 Heisenberg chain can serve as a perfect quantum 
channel. Furthermore, we apply the numerical method to mixed states and study the entanglement 
as a function of temperature. 

PACS numbers: 03.67.Mn, 03.65.Ud, 75.10.Pq, 73.43.Nq 
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I. INTRODUCTION 

The creation and distribution of entangled states plays 
a central role in quantum Information, because it is the 
key ingredient for performing certain quantum Informa- 
tion tasks, like teleportation 1] or quantum computation. 
In this respect, multiparticle quantum states can be con- 
sidered as entanglement resources, naturally appearing in 
many physical systems. On the other hand, it is believed 
that the study of multipartite entanglement might prove 
fruitful in other fields of physics, like condensed matter, 
e.g. for understanding the complex physics of strongly 
correlated states ^ In particular it has been shown 

lE lE lE Hl m that the ground state entanglement of 
various spin systems may exhibit characteristic features 
at a quantum phase transition. Hence it is desirable to 
find ways of characterizing and quantifying entanglement 
in multipartite systems. 

In 53 an entanglement measure, called the Localizable 
Entanglement (LE), has been defined. It quantifies the 
bipartite entanglement contained in a multipartite sys- 
tem. The concept of LE allows one to define the notion 
of entanglement length, which characterizes the typical 
distance up to which bipartite entanglement can be lo- 
calized in the System. Moreover, it has been shown for 
general pure qubit states that the LE can always be lower 
bounded by connected correlation functions jlOI] . Hence 
quantum phase transitions, characterized by a diverging 
correlation length, are equivalently detected by a diverg- 
ing entanglement length. In fact, the concept of LE has 
already proven useful in several studies of entanglement 
properties of spin systems 0, 0, 0, 0, Indeed, 
for ground states of many spin-1/2 models one finds that 
correlations and entanglement (as given by LE) typi- 
cally exhibit the same qualitative behavior. In gapped 



spin-1 systems, which exhibit finite correlation lengths 
[15j, however, a ground state has been found for which 
the entanglement length is infinite pdj . This example 
also shows that the concept of LE can serve to detect 
hidden order in certain states. 

The LE is the maximum entanglement that can be lo- 
calized, on average, between two parties of a multipartite 
System, by performing local measurements on the other 
parties. Hence, LE is defined in an operational way and 
has a clear physical meaning. For instance it can be 
used as a figure of merit to characterize the performance 
of quantum repeaters ^E|- Note that, in the context of 
LE, particles are not traced out but measured. This is 
in contrast to earlier approaches, where the concurrence 
of the reduced density operator of two separated spins 
in a spin chain has been calculated (see e.g. 0, IE|)- Al- 
though the concurrence exhibited characteristic features 
at a quantum phase transition, it does not detect long 
range quantum correlations. 

The fundamental difference between tracing and mea- 
suring can be illustrated with two simple examples: For 
both the GHZ state 17] and the cluster state ,1^ it can 
readily be checked that the reduced density operator of 
any two qubits contains no entanglement at ali. On the 
other hand one can find a local measurement basis, such 
that the LE is maximal. We further note that in the 
case of the GHZ-state the entanglement properties could 
also have been revealed by studying the connected ver- 
sion of the two-point correlation function. Indeed, the 
GHZ-state is one of many examples 10], for which corre- 
lations can be identified with quantum correlations, i.e. 
entanglement. However, this intimate conncction does 
not hold true for ali pure quantum states (as shown in 
[III]). For instance, in the case of the cluster state the 
LE is maximal, whereas the connected correlations are 
ali zero. 
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The intention of this paper is to provide a framework 
for using LE as an measure for localizable bipartite entan- 
glement of multipartite quantum systems. Starting from 
the original definition 10], we introduce several variants 
of LE, and establish basic rclations between them. Apart 
from the notion of entanglemcnt length , we also give a 
meaning to the notion of entanglement fluctuations in 
terms of LE. Moreover we generahze the earUer resuh, 
that the LE is lower bounded by connected two-particle 
correlation functions, to pure qutrit states, e.g. ground 
States of spin-1 systems. Apart from these analytical 
findings we present a method for the numerical computa- 
tion of the LE for arbitrary one dimensional spin systems. 
It allows US to efHciently simulate chains with even more 
than 100 sites and also works for finite temperatures. We 
apply this method for the study of LE present in ground 
states of various spins models. We find that the entangle- 
ment fluctuations as well as the LE exhibit characteristic 
features at a quantum phase transition. For instance, we 
observe a discontinuity in the first derivative of the LE 
at a Kosterlitz-Thouless transition |l9|. These findings 
can be understood as a direct consequence of the numer- 
ical observation, that for spin-1/2 systems the LE and 
the maximal connected correlation function are typically 
equal to each other. In the case of the spin-1 Heisenberg 
antiferromagnet, however, we rather observe the opposite 
behavior, namely an infinite entanglement length but fi- 
nite correlation length. In this context we also comment 
on a possible connection between the presence of hidden 
order in the state and the existence of long range en- 
tanglement. In particular, we present examples showing 
that, in general, such a connection does not exist. As 
a further application of our numerical method we study 
the effect of finite temperature on the LE. For the AKLT 
model [2Õj we find that the entanglement length increases 
exponentially with the inverse temperature. 

The paper is organized as foUows: We start by giv- 
ing a formal definition of LE and derive quantities like 
entanglement length and entanglement fiuctuations. We 
further connect the concept of LE with the idea of quan- 
tum repeaters. In Sect. III we recapitulate the earlier 
result for the lower bound of LE in terms of connected 
correlation functions and generalize it to pure qutrit sys- 
tems. Next, we present in Sect. IV a numerical method 
to compute the LE for spin chains. The numerical scheme 
is based on the matrix-product state (MPS) 0,112 rep- 
resentation of ground states and the Monte Carlo (MC) 
[2^ method. In Sect. V this method is applied to calcu- 
late the LE and the entanglement fluctuations for ground 
states of various standard spin-1/2 chains. In Sect. VI 
we study the LE of gapped spin-1 models. In particular, 
we present numerical calculations for the ground state of 
the Heisenberg antiferromagnet, showing that the entan- 
glement length diverges. Furthermore, we comment on 
a possible connection between a diverging entanglement 
length and hidden order in the system. In Sect. VII 
we demonstrate that our numerical method can also be 
used to calculate the LE of mixed states. As an example 



we compute the entanglement length as a function of the 
temperature for the AKLT model. In Appendix A we 
present an extended version of the proof in ^3 j leading 
to the lower bound of LE. We connect in Appendix B 
the entanglement of pure two qubit and two qutrit states 
to the maximum connected correlation function of those 
states. Finally, in Appendix C, we show how to calculate 
the LE (and the string order parameter) analytically for 
pure states, represented by MPS with qubit bonds. 



II. DEFINITION AND BASIC PROPERTIES OF 
LOCALIZABLE ENTANGLEMENT 

We consider a multipartite system composed of N par- 
ticles. With each particle we associate a finite dimen- 
sional Hilbert space. For simplicity we refer in the fol- 
lowing to the particles as spins. 



A. Definition 

The Localizable Entanglement (LE) of a multi-spin 
state p is defined as the maximal amount of entangle- 
ment that can be created (i.e. localized), on average, 
between two spins at positions i and j by performing 
local measurements on the other spins. More specifi- 
cally, every measurement Ai specifies a state ensemble 
'■= {PstPI-'}- Here ps denotes the probability to ob- 
tain the (normalized) two-spin state pl^ for the outcome 
{s} of the measurements on the N — 2 remaining spins. 
The average entanglement for a specific á4 is then given 
by: 



(1) 



where E{pl^) is the entanglement of pl^ . Suitable entan- 
glement measures will be discussed later in this section. 
The Localizable Entanglement is defined as the largest 
possible average entanglement: 



sup 
Mec ■ 



(2) 



with C denoting the class of allowed measurements. We 
call the measurement Ai which maximizes the average 
entanglement the optimal basis. It is important to note 
that the only restriction on Ai is that the measurements 
are performed locally i.e. on individual spins. Apart 
from that, the measurment basis is arbitrary and can also 
vary from site to site. We distinguish three classes C of 
measurements: projective von-Neumann measurements 
(PM), those corresponding to positive operator-valued 
measures (POVM), and general local measurements that 
allow also for classical communication of measurement 
results (LOCC). In terms of LE the foUowing relationship 
between these classes holds !24|: 



(p) < j Íp) < L-^. (p). 



(3) 
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In this paper we will be mainly concerned with projective 
measurements. To simplify thc notation wc omit in this 
case the superscript PM. 

The definition Q still leaves open the choicc of thc 
entanglement measure E for the states {pl^}- Suitablc 
measures depend on aspects hke the dimensionahty of 
the spins and the purity of the state. In the foUowing 
we specify for the cases of both pure and mixed spin-1/2 
and spin-1 systems appropriate entanglement measures 
used in this article. 

For pure multipartite states p, the states {pl^} after 
the measurements are also pure, and then there exists, 
in principie, a special entanglement measure |25|. i.e. 
the entropy of entanglement. For a pure bipartite state, 
Pij = \ip){'ip\, the entropy of entanglement Ee is defined 
as the von Neuman entropy of the reduced density oper- 
ator p, = trj(py ): 



EeÍPí]) = -tr(pi log2Pí). 



(4) 



In the case of a pure two qubit state, it can be shown 
that the Ee is a convex, monotonously increasing func- 
tion, Ee = f{C), of the concurrence [2^. The conur- 
rence C{iIj) for pure states and the convex function / are 
defined as: 



f{C) 



H 



1 + VI - 



H{x) := -xXog^x - {\ ~ x)\og2{l - x) 



(5) 
(6) 
(7) 



Here, denotes the complex conjugate of IV') in thc 
standard basis and H{x) is the Shannon entropy. In 
this article we will typically use the concurrence to mea- 
sure the entanglement of two qubits, because it can be 
simply related to connected correlation functions (see 
Appendix B). However, due to the convexity of the func- 
tion /(C), the LE as measured by the concurrence, L^, 
yields lower and upper bounds for the LE as measured 
by the entropy of entanglement, Lff: 



(8) 



Therefore the qualitative behavior of these two variants 
of LE will be very similar. 

In the case of mixed states the entanglement can be calcu- 
lated using, e.g., the entanglement of formation Ep H^l 
or the negativity |28j| . For two qubits an explicit formula 
for the Ep in terms of the concurrence C exists, which 
reduces to in the case of pure states. The negativity 
N can in principie be calculated for any spin dimension. 



B. Entanglement length and fluctuations 

In the field of strongly correlated systems and more 
specifically in the study of quantum phase transitions, 
the correlation length, is of great importance. The 



concept of LE readily lends itself to define the related 
entanglement length, Çe, as the typical length scale at 
which it is possible to create Bell states by doing local 
measurements on the other spins: 



lim 

n^oo \ n 



(9) 



The entanglement length is finite iff the LE 
Lfi^n exp(— ri/^s) for n 00, and the entan- 
glement length ^E is defined as the constant in the 
exponent in the limit of an infinite system (see also 
Aharonov 

Let US now have a closer look at the statistical nature 
of LE, as it is defined as an average over ali possible mea- 
surement outcomes (see For practical purposes one 
can only control the measurement basis but not a spe- 
cific outcome. Therefore it would be useful to have an 
estimate of how much the entanglement of a particular 
measurement outcome deviates from the mean value as 
given by the LE. This Information is contained in the 
variance of the entanglement remaining after measure- 
ments. We can thus define the notion of entanglement 
fluctuations: 



{5L 



M,E\2 



(10) 



The entanglement fluctuations can be defined for any 
measurement M. Typically we choose for M the op- 
timal basis, which maximizes the average entanglement. 
In this case we drop thc index M. in UlUII . 

The study of both the entanglement length and the 
entanglement fluctuations could provide further inside 
in the complex physics of quantum phase transitions by 
revealing characteristic features at the quantum criticai 
point. Examples for this are presented in Sect. V. 



C. Connection to quantum repeaters 

So far we have given a purely mathematical definition 
of LE (|2Jl. However, it is evident that the LE is defined 
in an operational way that can directly be implemented 
on certain physical systems. In addition the concept of 
LE may also play an essential role in some interesting 
applications of quantum Information theory. To be more 
precise, LE can serve as a figure of merit for the "perfor- 
mance" of certain kinds of quantum repeaters (QR). 

Many tasks in quantum Information processing require 
long-distance quantum communication. This mcans 
quantum states have to be transmitted with high com- 
munication fidelity via a quantum channel between two 
distant parties, Alice and Bob. Since quantum transport 
is also possible via teleportation this problem is equiv- 
alent to establishing nearly perfect entanglement between 
two distant nodes. AU realistic schemes for quantum 
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communication are presently based on the use of pho- 
tonic channels. However, the degree of entanglement gen- 
erated between distant sites typically decreases exponen- 
tially with the length of the connecting physical channel, 
due to hght absorption and other channel noise. To over- 
come this obstacle the concept of quantum repeaters has 
been introduced The central idea is to divide the 

channel into segments and to include additional nodes. 
Entanglement between adjacent nodes can be extended 
to larger distances using entanglement swapping foUowed 
by purification. After several rounds one obtains a pair of 
almost maximally entangled nodes, shared by Alice and 
Bob, that can be used for perfect quantum transport via 
teleportation. A possible physical realization of the QR 
using trapped atoms is sketched in Fig. ITl [8?ÍI|. 



Node 1 



Node2 



B 




not have an explicit parameterization of the state un- 
der interest, but just Information about the classical one- 
and two-particle correlation functions (which allows one 
to parameterize completely the two-spin reduced density 
operator pij). It would therefore be interesting to derive 
tight upper and lower bounds to the LE solely based on 
this Information. 



A. Upper bound 

The upper bound can readily be obtained using the 
concept of entanglement of assistance (EoA) (3l]|. The 
EoA depends only on the reduced density operator pij 
for spins i and j and can be defined for any spin dimen- 
sion. In the case of a pair of qubits an explicit formula 
for EoA can be derived 1^3 • Given pij and a square 
root X, Pij = XX'^ , then the EoA as measured by the 
concurrence reads: 



i?i^"(p.,):=tr|X^K®aJ)X| 



(11) 



FIG. 1: Illustration of the quantum repeater scheme for 
trapped atoms connected by optical fibres. Qubits are rep- 
resented by the internai states of the atoms. Applying laser 
beams the internai states of atoms in adjacent cavities become 
entangled via the transmission of photonic states. CoUective 
measurements on the nodes (indicated by arrows) foUowed 
by purification lead to the generation of a nearly perfectly 
entangled pair of qubits between A and B. 

Let US now discuss how a QR setup can be character- 
ized by the LE. First of ali, let us identify the particles 
sitting at different nodes by spins. It is important to note 
that, by combining several spins to a larger Hilbert space 
of dimension d, the operations required for purification 
and entanglement swapping on this set of spins can be in- 
terpreted as local operations on a single spin of dimension 
d. Thus the QR can be treated as a system of interacting 
spins being in a state p. In order to assess and quantify 
the usefulness of such a setup as a QR one has to com- 
pute the foUowing figure of merit: What is the maximum 
amount of entanglement that can be generated between 
the two end spins by performing local operations on the 
intermediate spins? But this number is nothing else than 
the LE. The question, which variant of LE (j2Jl should be 
used, depends on the class of available local operations 
(PM, POVM or LOCC). Typically classical communica- 
tion is allowed so that L^^^*-^ has to be taken as figure of 
merit. However, not every measurement might be physi- 
cally realizable. Therefore the LE will in general give an 
upper bound for the performance of a given QR setup. 



III. BOUNDS ON LE 

Due to its variational definition, the LE is very diíR- 
cult to calculate in general. Moreover, typically one does 



with \A\ 



where 
i^^X. 



' (7: 



Viu. Hence (p.,,) = ELi<^k, 
Gk are the singular values of the matrix AT (cr* -j^ 
Note that for pure four qubit states a variant of EoA with 
local measurements was considered in |32j. 



B. Lower bound 

First, from the definition of LE ^ it foUows that any 
specific measurement, e.g. in the computational basis, 
trivially provides a lower bound on the LE. More inter- 
estingly, it has been proveu in [TÕI | that for general pure 
qubit states the LE can be lower bounded by connected 
correlation functions. In the foUowing we recapitulate 
the central results leading to this bound for spin-1/2 sys- 
tems and present some extensions. Next we show to what 
extent these findings can be generalized to higher dimen- 
sional spin Systems. 



1. Spin-1/2 Systems 

The basic idea is to establish a connection between the 
LE and connected correlation functions of the form: 



Q 



AB 



tr[p( S\ ® Si 



tr[p( S\ ® l)]tr[/a(l - 



.S^)]. (12) 



For qubits the operators Sa, Sb can be parameterized 
by directions a, b, representing unit vectors in a 3D real 
space: Sa = a ■ a, Sb = b ■ a with a = {ux, cr„ , Q-y ) . 
We start out by quoting the central result of |10| : 
(l.i) Given a (pure or mixed) state of N qubits with con- 
nected correlation function Q^^g between the spins i and 
j and directions a, b, then there always exists a basis in 
which one can locally measure the other spins such that 
this correlation does not decrease, on average. 
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An extended version of the proof is presented in Ap- 
pendix A. There we also show that this result can be 
generalized to a setup, where the spins i and j can be 
of any dimension, but the remaining spins (on which the 
measurements are performed) are still qubits. In a spin- 
1/2 System such a situation can arise, for example, when 
considering correlations between two blocks of spins. 

Next, we relate correlations with entanglement. We 
note that after the measurement process and for an ini- 
tially pure state we end up with a pure state of two 
qubits. For such a state we have proven the foUowing 
result ^3 (see also Appendix B): 

(l.ii) The entanglement of a pure two qubit state \ipij) 
as measured by the concurrence is equal to the maximal 
correlation function: 

C{^,,)=mj,x\Q'Í^{^,j)\. (13) 

ã.b 

Combining (l.i) and (l.ii) we know that for a given pure 
multi qubit state \^p) and directions a, b there always ex- 
ists a measurement M such that: 

QÍbW < T.P^QabO^s) < Y^PsCm- (14) 

s s 

The term on the very right is equal to the average entan- 
glement as measured by the concurrence L^'*^, which 
trivially is a lower bound to the LE as defined by L^. 

Since the directions a, b can be chosen arbitrarily, rela- 
tion H14|l holds in particular for directions maximizing 
*3as (^)- Hence we can establish the desired lower bound 
on LE li^: 

(l.iii) Given a pure state \il>) of N qubits, then the LE as 
measured by the concurrence is larger or equal than the 
maximal correlation: 

Lf^.(V;)>max|Q:^'5(^)|. (15) 

Making use of the basic properties of LE, presented in the 
previous section, we can immediately derive analogous 
bounds for some other variants of LE; for example, 

if;(^) > fimB^^m^m). (16) 

^POVM,E.(^) > /(m_ax|Q:^"^(^)|), (17) 

with / being the convex function defined in J^. We will 
see below that relation (|17|l can be generalized to spin-1 
Systems. 

2. Higher dimensional spin systems 

We now try to extend the previous findings beyond 
spin-1/2 systems. First, we look for a generalized version 
of statement (l.i). Unfortunately the techniques used in 
the proof for qubits seem to fail already for qutrits. Nev- 
ertheless, a generalization is still possible by changing a 



little bit the perspective. For this we embed a spin-S* 
in a higher dimensional Hilbert space, being composed 
of n > log2(2iS' + 1) virtual qubits. Let us denote the 
(25'+ 1) X 2" matrix governing this transformation by P. 
In the case 25* + 1 = 2" the embedding is trivial and 
the situation becomes equivalcnt to the qubit case. Thus 
the result (l.i) can immediately be generalized, because 
local measurements on the virtual qubit systems can be 
chosen such that (l.i) holds. 

In the case 25 + 1 < 2" a similiar argument applies if 
we allow for POVM measurements on the spin- 5 sys- 
tem. To be more precise, let us consider a mixed state 
p of three spin-S* particles. The spin on which the 
measurement is performed (let us denote it with the 
Índex 3) is embedded in a 2" dimensional system. The 
embedded state is then given by the transformation: 

p' = (I12 «1 P^)p{li2 «1 P3). In the Hilbert space of the 
n virtual qubits one always finds local projective mea- 
surements {Maj^...a„} = {\oii){ai\ (8) ... (8) |Qí„)(a„|} such 
that (the generalized version of) (l.i) holds for the state 
p' . In terms of the original state p this measurement in 
the 2" dimensional space corresponds to a POVM mea- 
surement {PMa^ a„P^} on a spin-S system, because: 

Ea,...a„ PMo.,...o.J^ = P12"X2"P^ = 1(2S+1) X (2S+1) • 

Thus one can generalize the result (l.i) to arbitrary 
spin dimensions in the foUowing way: 
(2.i) Given an arbitrary multi-spin state with connected 
correlation function between spins i and j for ar- 

bitrary operators Sa^Sb, then there always exists a local 
POVM measurement on the other spins such that this 
correlation does not decrease, on average. 
We note that the lower bound in (2.i) can already be 
reached by applying local measurements on the virtual 
qubit system. Performing joint measurements (e.g. Bell 
measurements on pairs of qubits as shown in ,11]) can 
lead to a considerable enhancement of the average corre- 
lations. If in addition 2" — (25 -f 1) joint measurements 
can be chosen such that they are orthogonal to the pro- 
jector P, the resulting measurement on the spin-S system 
corresponds to a projective von-Neumann measurement. 

The most difficult part is to establish a connection be- 
tween correlations and entanglement for pure two spin 
States in analogy of (l.ii). In the case of qubits we made 
explicitly use of the fact that the group SU{2) is the cov- 
ering group of 50(3). Moreover the concurrence served 
as an entanglement measure, which was easy to handle. 
For higher spin dimensions we refer to the entropy of en- 
tanglement Ee 101) as a suitable entanglement measure 
for pure bipartite states. In the special case of qutrits we 
were able to show the foUowing relation (Appendix B): 
(2.ii) The entanglement of a pure two qutrit state 
as measured by the entropy of entanglement can be lower 
bounded by: 

i?ií(V'.,)>/(max| Q%[^^j)\), (18) 

where / is the convex function © and Sa, Sb in Q^ab 
(|12|l are operators, whose eigenvalues lie in the interval 
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[— 1; 1]. Combining again (2.i) and (2.ii) we can formu- 
late a bound on LE for spin-1 systems: 
(2.iii) Given a pure state of N qutrits, then the LE 
as measured by the entropy of entanglement and which 
allows for POVM's, is lower bounded by the maximum 
connected correlation function in thc foUowing way: 

LfO^^'^-(^)>/(max|QÍ^"5(V^)|). (19) 

In summary, we have shown for pure qubit and qutrit 
states that connected correlation functions provide a 
lower bound on LE. This bound allows for two intrigu- 
ing limiting case: (i) Entanglement and correlations may 
exhibit similiar behavior. (ii) Spins may be maximally 
entangled although they are uncorrelated in the classi- 
cal sense. In thc forthcoming sections we will present 
examples for both scenarios. 



Without loss of generality we can assume that the com- 
putational basis is the optimal one. Otherwise we can 
make a change of basis in H21|) . Hence the LE is given 
by: 

L.^(V'mp)=5]pm£;(0m) , (22) 

{s} 

where p[s} = ((/«{s} |'/'{s})/(^mp|V'mp) is the prob- 
ability for obtaining the normalized state l^is}) — 

l'/'{s})/('/'{s}l0{s})^^^- Though the MPS representation 
allows US to efficiently calculate the states after the mea- 
surements we still face the problem that the sum (|22|) 
involves an exponential number of terms {d^~'^). To find 
a good approximation of this sum we propose a scheme 
based on the Monte Carlo (MC) method, which will be 
now explained in more detail. 



IV. COMPUTATION OF LE BASED ON 
MATRIX PRODUCT STATES AND THE MONTE 
CARLO METHOD 

In this section we propose different techniques to cal- 
culate the LE numerically for ground states of ID spin 
systems. Lct us consider a chain of N spins of dimen- 
sion d — 2S + 1. The ground state of the system can 
be determined exactly by diagonalization of the Hamil- 
tonian. However, the state is characterized by an expo- 
nential amount of parameters as a function of A'', thus 
limiting the exact treatment of the problem to relatively 
small system sizes. 

As an alternative approach we can start out with an 
approximation of the exact ground state in terms of the 
the so-called matrix product states (MPS) (21i.i22l]: 

d 

|^MP)= Tr(Ar ...A^")|si,...,SAr). (20) 

si ,...,sjv=l 

Here the state is described by N matrices A^^ of dimen- 
sion D. We note that the MPS is written in the 
computational basis and accounts for peri odic boundary 
conditions (PBC). It has been shown |2l lã^ that MPS 
appear naturally in the context of the density-matrix 
renormalization group (DMRG) method "s^. Assuming 
we are able to calculate the ground state in MPS form 
(|20|l , let US now present a scheme to compute the LE from 
that. For translationally invariant systems it is sufíicient 
to consider the LE between the spins 1 and j = í + n. 
The (pure) normalized state of these two spins, after lo- 
cal measurements on the remaining ones have been per- 
formed, is conditioned on the measurement outcomes de- 
noted by the (N-2)-tuple {s} := {s2 ■ ■ . Sj-iSj+i . . . sn} 
and proportional to 

d 

= (MI^MP) = E Tr(Ar...A?,«)|si)|s,). 

Si,Sj = l 

(21) 



A. Monte Carlo method 

The Monte Carlo method provides an efficient way 
of selecting M states sequentially from the (given) 
probability distribution {p{s}}- Thc LE can thus be ap- 
proximated by: 

1 *^ 1 

^fMC),y- (V'MP)«^ 5] i?(0^) ± -=SLf<iM. 

(23) 

Note that the accuracy of the MC method depends on 
the entanglement fluctuations SLfj, which can be com- 
puted within the MC scheme as well. 
For selecting the states we foUow the Metropolis 
algorithm (23i | and use single-spin-flip dynamics. We 
start with an initial state j^^) corresponding to a spe- 
cific measurement outcome {s}. /From this we create a 
trial state {(jji,) by randomly picking a site i and changing 
the state of this spin with equal probability according to 
Si — > mod(si ± 1, d) for Si = 0, 1, ... d — 1. For a spin-1/2 
this simply amounts to a spin flip. The probability p^, for 
obtaining the trial state after a measurement can conve- 
niently be computed using the MPS representations l|21(l . 
The trial state is accepted with probability: 

P(.^.).{f (24, 

If thc trial state is accepted it serves as a starting point 
for creating a new trial state. After N — 2 steps, defining 
one MC sweep, the entanglement of the current state is 
calculated |2a|- After M sweeps the algorithm stops and 
the average (|^ is performed. 

B. Finding the optimal measurement basis 

The definition of the LE (|2Jl requires an optimization 
over ali possible measurement strategies. A good guess 
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for the optimal basis can typically by found using exact 
diagonalization for small system sizes, foUowed by numer- 
ical maximization of the average entanglement. Alterna- 
tively, the optimal basis can also be extracted directly 
from the MPS matrices Ai in H21() . Using a generahzation 
of the concurrence for pure bipartite D x D states, it has 
been shown in (llj (for an open chain with ZJ-diniensional 
spins at the ends), that the optimal basis is the same 
basis, that maximizes the expression 'J2s- I det(A^')p/^. 
Hence we can consider the matrices Ai as (unnormalized) 
pure D X D xd states, for which we want to calculate the 
LE with respect to the D x D system. Since we measure 
only on a single site this problem is equivalent to calcu- 
lating the EoA of the reduced (1?^ x D^) density matrix, 
which can be done numerically (see also Appendix C). 
For the models we studied our numerical analysis indi- 
cates that these findings hold independent of the choice of 
both the entanglement measure and the boundary condi- 
tions. As a further numerical result we find that the op- 
timal basis appears to be independent of the system size. 
Hence exact diagonalization and numerical optimization 
for small systems usually provides the most efficient way 
to find the optimal measurement strategy. 



C. Determination of the MPS 

In the foUowing we are interested in finding a good ap- 
proximation of the true ground state in terms of MPS 
(|20|l for systems with PBC. We choose PBC in order to 
minimize boundary effects and to better mimic the be- 
havior in the thermodynamic limit alrcady for small sys- 
tem sizes. Let us now introduce a method to determine 
translationally invariant MPS. It is based on the follow- 
ing idea: We use the DMRG algorithm for an infinite 
chain to extract a site- independent set of matrices A*, 
defining a translationally invariant MPS l|2(J|) for infinite 
N. We then us the same set A" to construct a MPS with 
PBC for arbitrary A''. It is obvious that this method is 
very efhcient, particularly for large A^^, because we have 
to run the DMRG only once to obtain the MPS repre- 
sentation for any system size. One might expect that 
this increase in efficiency happens at the cost of preci- 
sion. However, an optimum in accuracy on the part of 
the MPS is not crucial for the numerical calculation of 
the LE, since the limiting factor for the accuracy is typ- 
ically the MC method. 

We start by briefly reviewing the variant of DMRG, 
represented hy B • B [s^, for an infinite ID chain. At 
some particular step the chain is split into two blocks and 
one spin in between. The left block (i) contains spins 
1, . . . , M - 1, and the right one (R) spins M + 1, . . . ,N. 
Then a set oi D x D matrices A"'^" is determined such 
that the state 

d D 

I*) = E E ® \s}m ® (25) 

s=l a, (3=1 



minimizes the energy. The states \a)L.R are orthonormal, 
and have been obtained in previous steps. They can be 
constructed using the recurrence relations 

D d 
q' = 1 s=1 

where the block L' contains the spins 1, . . . , M — 2. 

Numerically we find that the matrices {/['^l''* can be 
chosen in such a way that they converge to (site inde- 
pendent) matrices at the fix point of the DMRG algo- 
rithm. Applying an appropriate transformation R, these 
matrices can be used to construct a translatio nally 
invariant MPS ^ with A' = RWR^^ = RÃ^R^ 0. 
Note that starting from these states, it is possible to cal- 
culate expectation values of products of local observables 
miillHi, since 

(vI/|Oi...O^|vI/)=.Tr(i?W (27) 

where 

d 

^o'= E {s\Ok W)Al^(AÍy . (28) 

s,s' — 1 

We applied this method to compute ground state expec- 
tation values using expression (|27|l for various ID spin 
systems with PBC f3q|. Our numerical analysis of sys- 
tems with finite N shows that both the energy and the 
correlations can be computed rather accurately. The 
achieved accuracy is several orders of magnitude higher 
than finite size effects, but also several orders of mag- 
nitude lower compared to the variational method [35|, 
which has been introduced recently. These findings even 
hold for system sizes as low a A « 10, and are rather 
surprising, because the MPS is constructed from the in- 
finite chain. We further checked that also the long range 
behavior of correlations and entanglement is reproduced 
correctly by our translationally invariant MPS. 

To sum up, our numerical results indicate that transla- 
tionally invariant states can be sufficiently well approx- 
imated by a single set of MPS matrices A^ for almost 
arbitrary system size A^. This observation, together with 
the findings in [TH . might shed some light on our pre- 
vious numerical finding that the optimal measurement 
basis for LE is typically both site-and size- independent. 



V. EXAMPLES OF LE IN SPIN-1/2 MODELS 

In this section we apply the concept of LE to quan- 
tify the localizable ground state entanglement of various 
spin-1/2 models. After some general considerations we 
compute the LE as measured by the concurrence Lfj nu- 
merically for two specific examples. 
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A. General considerations 

We consider spin-1/2 Hamiltonians of thc form 

H=-Y. E 1>>Í-Y.^'^^' (29) 

i.j a—x.y,z i 

with parity symmetry, [ií, H^] — O and := (Xi^icr^- 
Extcnsive numerical calculations on systems of up to 20 
qubits showed that our lower bound is always close to 
the LE as measured by concurrence Lfp and typically is 
exactly equal to it: this is surprising and highlights the 
thightness of the given lower bound. Note also that when- 
ever parity symmetry is present, the upper and lower 
bound are given as foUows [TÕ| : 

max(|Q^J,|Q^J,|g;i|) < Lf^^ < y^±—l— , (30) 
4 = (l±(aX))' - {{al)±{ai)f. 

The fact, that the lower bound is usually tight, can 
also be derived from the numerical observation that for 
Hamiltonians of the form (|29|l measurements in the (stan- 
dard) CTz-basis {A4 = Z) yield in most cases the opti- 
mal result. Expanding the ground state in that basis, 
\ip)o = X] 1*1 • ■ ■''■n), it is straightforward to show 
that e.g. in the case of constant (site independent) cou- 
pling and 7j, the ground state energy is minimizcd 
if ali expansion coefficients Cíj...í„ have the same sign. 
This guarantees, together with the parity symmetry, that 
the average entanglement Lfj~^ for measurements in the 
standard basis is equal to either the x x or y — y cor- 
relation. To be more precise we distinguish the foUowing 
cases: 

(7^~7^)(7^+7^)>0: hff ^\{alai)l (31a) 
(7^-7-')(7^+7^)<0: hff ^ \{alal)l (31b) 

where the conditions refer to ali sites i and j of the 
chain. Most of the prominent spin Hamiltonians studied 
in literature, like the Heisenberg, XY or XXZ model etc, 
trivially fulfill one of the conditions (|31|l . because their 
coupling coefficients are site-independent. Hence, mea- 
surements in the standard basis would yield localizable 
quantum correlations that are completely determined by 
classical correlations. 



B. Ising model 

As an illustration, let us now discuss the LE 
of the Ising model in a transverse magnetic field 
{la = A ^f> r t^i.i+i ; 7' = 1 in (P^ ). which has been solved 
exactly j39j and exhibits a quantum phase transition at 
A = 1. In this case, the maximal connected correlation 
function is always given by Qxx^ which thus yields the 
best lower bound on LE. Numerical optimization for a fi- 
nite chain indicates that the standard basis is indeed the 



0£timal one and thus the lower bound is equal to Lp^- 
|lOi | . We checked analytically, using perturbation theory, 
that for an infinite chain this numerical result is indeed 
true. However, for a spin distance n — \i — i\ one has to 
go to n-th order perturbation theory, limiting this ana- 
lytical treatment to rather small n. 

Consequently, we can use exact results for the connected 
correlation function Q^x ^âj to completely characterize 
the behavior of the LE in the Ising chain. The Ising 
system is therefore also an ideal candidate for testing 
the performance of our numerical method outlined in 
Sect. IV. In Fig. Hwe plot L,^,+„ and Q^^*+" as a func- 
tion of the spin distance n for a chain with iV = 80 sites. 
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FIG. 2: Calculation of the I/pi+,i(diamonds) for the 
ground state of the Ising chain in a transverse field 
( 7a = A5cí,a;5j,i+i; 7* = 1 in II29II ') as a function of the spin 
distance n. For comparison we plot the exact result [39l| for 
the correlation function Q^'^^"'. Left: A = 0.8, exponential 
decrease; Right: criticai point A = 1, power law decrease 
(~ n~^^*); Numerical parameters (see Sect. IV): A'^ — 80, 
D = 16, MC sweeps M = 20, 000. 

For A < 1, the LE decreases exponentially with n, and 
the entanglement length is finite. At the quantum criticai 
point A = 1, the behavior of the LE changes drastically, 
because it suddenly decreases as a power law, Lf^^^ ~ 

n~^/'^, thus leading to a diverging entanglement length 
^B. In Fig. [21we observe that the MC method becomes 
less accurate at the criticai point. As we will see later, one 
reason is that the statistical error due to entanglement 
fluctuations becomes rather large at the criticai point (see 
Fig. O. Another (systematic) error might be induced 
by the single-spin-flip dynamics used to create the trial 
state. Better results for the criticai region could possibly 
be achieved by applying the Wolff algorithm [41|. Here, 
a cluster of spins depending on their spin orientation is 
flipped, which accounts for the formation of domains. 
In Ji)J it was shown that for the case A > 1 we also get 
= oo, since the LE saturates to a finite value given 
by Ml = 1/4(1 - A-2)i/4, Indeed, the ground state is 
then close to the GHZ-statc. In a more realistic setup, 
however, the parity symmetry of the Ising Hamiltonian 
will be broken by a perturbation and the ground state 
for large coupling will also be separable, as it is given by 
a superposition of two GHZ states with different parity 
0. 
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Let US now study the behavior of entanglement fluc- 
tuations <5L£- inherent to the statistical definition of LE 
(see Def. (|10|l l as a function of the couphng A. In Fig. El 
we plot SLfj for different parameters n and N. 




X X 



FIG. 3: Exact calculation of entanglement fluctuations 
(5-í/i^í+„ liou as a function of the couphng parameter A for 
a finite Ising chain in a transverse magnetic held with PBC. 
Left: distance n = 1, iV = 6 (dashed), iV = 12 (dotted), 
N ^ 16 (solid); Right: distance n = 4, N = 9 (dashed), 
iV = 12 (dotted), iV = 16 (solid). 

The maximum of the fluctuations is always located in 
the vicinity of the criticai point A = 1 and gets shifted to 
larger A values with increasing N. Thus the increasing 
entanglement fluctuations reflect very well the increas- 
ing complexity of the wavefunction close to the criticai 
region. The location of the maximum A™ in the thermo- 
dynamic limit (iV —^ oo) apparently depends on the dis- 
tance n of the two spins. For nearest neighbors (n = 1) 
we observe that the maximum of õLf^^-^ is somewhat 
shifted to the right of the criticai point (Am ~ 1.025). 
For ali distances n > 1, however, our numerical calcula- 
tions show that the maximum is positioned at A < 1 but 
becomes asymptotically close to the criticai point with in- 
creasing n (for n — A see Fig. O. Furthermore, in Fig. O 
we see that the absolute value of the maximum increases 
with n and becomes comparable with Lfj itself (42] . The 
strong fluctuations inherent to the Ising model lead to 
large statistical errors in the numerical calculation of LE 
using Monte Carlo (see (I23|l 'l. The errors become even 
more pronounced for the calculation of the fluctuations. 
This is the main reason, why we have restricted ourselves 
here to exact calculations for a small system with PBC. 
However, we confirmed that the data for A'^ = 16 repre- 
sents the behavior in the large N limit reasonably well 
and no qualitative changes occur. 



C. XXZ model 

Let US now turn to the discussion of another ex actlv 
solvable ID spin system, the so called XXZ model |43) . 
This model not only appears in condensed matter physics 
in the context of ferro- or antiferromagnetic materiais. 
Recently it has been shown that it can also effectively 
describe the physics of ultra cold atoms in a deep optical 



lattice |43l. The Hamiltonian can be written as 

HxKZ = - cr^"^^ + ^1 (^'y^ + A o-^ 0-^+^ -I- jcr^ ], (32) 

i 

where we have introduced two dimensionless parameters, 
which can be varied independently: the anisotropy A and 
the magnetic field h/ J va units of the exchange coupling. 
The phase diagram of the XXZ model as a function of 
these two parameters is depicted in Fig. 0] 




FIG. 4: Schematic drawing of the phase diagram of the XXZ 
model (Eq. Il32ll 'l as a function of the anisotropy A and the 
magnetic field h/ J |4Í^ |. In regions A and C the ground state 
has an energy gap, whereas in region B the system becomes 
gapless (criticai) . Point E is the ferromagnetic XXX point and 
point D corresponds to the antiferromagnetic XXX point. 

The XXZ model can be solved exactly using the Bethe 
ansatz J^J. Unfortunately analytical expressions for the 
correlation functions, which would yield lower bounds for 
LE, have only been worked out in special cases. E.g. 
for the antiferromagnetic XXX model in a magnetic field 
(line of constant A = — 1 in Fig. Q arialytical solutions 
for the correlations are summarized in jl2l | . Numerics on 
a flnite chain of up to 14 spins show that again measure- 
ments in the cr^-basis appear to be optimal. Using the 
result (|31a|l this implies that L'(íj_^ — {a^a^"^). Hence 
the bounds given in [l^ are tight. In particular this 
means that the entanglement length is zero for h/ J > A, 
and infinite for h/ J < 4. 

After this introductory remarks let us study the ground 
state entanglement properties of the XXZ model in more 
detail. In particular, we are interested in finding charac- 
teristic features in the LE at the quantum phase transi- 
tions indicated by the lines a and ò in the phasediagram 
Fig. ^ For this purpose we calculate in the foUowing the 
LE (for flxed n) numerically, using exact diagonalization, 
as a function of the two parameters A and h/ J. 

1. LE as a function of the magnetic field h/ J 

In Fig. O we plot and the lower bounds Q^^^ 

and Q"^^^ as a function of the fleld h/ J for flxed A — 0.5. 



10 



We find that i^í+„ = (o-^cr*+"). This result can be un- 
derstood as a consequence of the numerical observation, 
that the standard basis appears to be optimal in the en- 
tire half space A > — 1. At the criticai point the LE 

1.2| ' ■ ' 1 

0.8 ^ "\ 




h/J 



ter suitable projective measurements in the z-direction 
these states reduce to maximahy entangled Beh 
States lí-^) = l/%/2 (|01) ± |10)). Hence in this 
limit Lfj — Ql^^ — ^ 1 and the cra;-basis turns out to be the 
optimal one. This line of reasoning strictly holds only in 
the limit A — oo, but can qualitatively be extended to 
the whole region C (Fig. 0)). 

For h/J > 1.6 we observe in Fig. [3 one of the rare 
examples, for which the LE is not exactly equal to the 
maximum correlation function. Note also that in this re- 
gion the maximum correlation function changes from Q!/^ 
to When approaching the criticai point h/J kí 6 

the standard basis becomes close to optimal again and 

In the foUowing we have a closer look on the phase tran- 
sition between the regions B and C in Fig. ^ 



1 
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2. LE as a function of the anisotropy A 

We consider the case of zero magnetic field and study 
the LE as a function of A in the vicinity of the antiferro- 
magnetic XXX point. At the criticai point A = — 1 the 
ground state undergoes a Kosterlitz-Thouless quan- 
tum phase transition. In Fig. Elwe have calculated nu- 
merically the localizable entanglement Lf^^^ for nearest 
neighbors and the corresponding lower bounds Q^'^"*"^ and 



FIG. 5: Calculation of the localizable entanglement -Í/í^í+„ 
(diamonds), the upper bound given by the EoA (IIH (dashed) 
and the lower bounds Q":^" (solid) and Q"^" (dotted) as a 
function of the field h/J for the ground state of the XXZ 
model Il,'i2|l . The numerical calculation is performed using 
exact diagonalization of a chain with N = 16 sites and PBC. 
The distance of the two spins is n = 4 sites; Up: A — 0.5; 
Down: A = -2; 



becomes zero, because the ground state is given by a 
product state. Note that this phase transition (indicated 
by line a in Fig. is sharp even for finite systcms, since 
it is due to levei crossing |43 ■ 

We now investigate the region A < — 1 , which contains 
a second quantum phase transition (indicated by line b 
in Fig. 21). The dependence of the LE and its bounds on 
the magnetic field for fixed A — —2 is depicted in Fig. O 
We again observe a sharp phase transition at h/J « 6, 
when the system enters the unentangled phase A (Fig. 

At /l/J «1.6 the LE experiences a sudden drop- 
off. However, it remains to be checked, whether this step 
is due to finite size effects or a characteristic feature of 
a quantum phase transition. Note that for h/J < 1.6 
we obtain that Lf^ = Q^J^- This feature can be un- 
derstood by considering the limiting case of zero field 
and A — > — oo. There, the Hamiltonian iíxxz commutes 
with the parity operator in x-direction, II^ — ^iLi'^x^ 
and the (doubly degenerate) ground state for even N 
is given by jV') = 1/V2 (|0101...01) ± |1010...10)). Af- 
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FIG. 6: Calculation of Lfi^i (diamonds) and the lower 
bounds Q"x^ (solid) and Q"^^ (dotted) as a function of the 
anisotropy A for the XXZ model <\'A2^ with zero field h/J. 
At the criticai point A = — 1 the LE exhibits a cusp, in con- 
trast to the monotonic behavior of the correlation functions. 
The calculation is performed using exact diagonalization and 
numerical optimization for a chain of length TV = 10 with 
periodic boundary conditions. 



One sees that Lf^i+i 



is equal to the maximum corre- 
lation function. However, at the criticai point the max- 
imum correlation function changes, due to the crossing 
of Q^'^^^ and Qyi'^^, thus leading to a cusp in the LE. 
Hence the quantum phase transition is characterized by 
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a discontinuity in the first derivative of ip^+i- This re- 
sult is remarkable, because for this Kosterlitz-Thouless 
transtion the ground state energy (i.e nearest neighbor 
correlation functions) and ali of its derivatives are conti- 
nous |45j . As shown in [ÍgI li^ , the concurrence (or the 
derivative) of the reduced density matrix pi^i+i also does 
not exhibit a discontinuity at the criticai point. This 
result is expected, because the concurrence is a function 
of one-particle and two-particle correlation functions. We 
further note that, according to our numerical analysis the 
reason for the cusp in is that the optimal measure- 

ment basis changes at the criticai point abruptly from the 
Gz- to the -basis. 

So far this is a purely numerical finding for a finite 
dimensional system. However, given that the average 
entanglement can be maximized by applying the same 
unitary transformation on ali spins, one can rigorously 
show that indeed a cusp in the LE must occur exactly 
at the criticai point and independently of the size N . 
The argument goes at foUows: At the antiferromagnetic 
XXX point the Hamiltonian (|32|l possesses SU{2) sym- 
metry. This means that any measurement basis yields 
the same LE. In particular, we know from (|31a|l that for 
measurements in the standard basis the LE is equal to the 
correlation function Q^'^"'"^. At the criticai point we thus 
have: ifí+i = Qzl'^^ = QxV'^- Since connected corre- 
lation functions yield a lower bound to Lfi+i (see Sect. 
III) the localizable entanglement Lf^_^_^ must exhibit a 
cusp at the criticai point, where Q^']^^ and Q^'^"*"^ cross. 

In summary our discussion of the LE in various spin- 
1/2 models has shown in which parameter regimes these 
systems can be used for e.g. localizing long-range en- 
tanglement as indicated by the entanglement length. We 
have further seen that the study of LE as well as the 
entanglement fluctuation provides a valuable tool for de- 
tecting and characterizing quantum phase transitions. In 
addition our numerical results indicate that the ground 
state entanglement of spin-1/2 Hamiltonians with two- 
spin nearest neighbor interactions is typically very well 
described by the maximum correlation function. We will 
now see that this observation is not necessarily true for 
spin-1 systems. 



VI. LE IN THE SPIN-1 
ANTIFERROMAGNETIC HEISENBERG CHAIN 

We study the ground state entanglement of the (gen- 
eralized) antiferromagnetic Heisenberg chain: 



Hai 



N~l 

El 



Si ■ Si- 



(33) 



which includes a biquadratic term. In a recent work [5Õ| 
it has been demonstrated that quantum Hamiltonians 
of this kind can be implemented with ultra cold atoms 
trapped in an optical lattice potential. 



A. General considerations 

For the Heisenberg antiferromagnet (AF) {f3 — 0) it 
has been shown by Haldane [sj that in the case of half- 
integer spins the spectrum is gapless in the thermody- 
namic liniit, and thus the correlation length of the ground 
state is infinite. For integer spins, however, an energy gap 
emerges, resulting in a finite correlation length. Let us 
now investigate whether a similiar connection holds for 
the LE and the corresponding entanglement length. 
^From the lower bound H15|l it foUows that the predicted 
infinite correlation length in spin-1/2 systems (with inte- 
ger log2(25'-|-l)) automatically implies a diverging entan- 
glement length. However, in the case of integer spins the 
correlation length is finite. Hence the lower bound 1)16(1 
for spin-1 systems includes the intriguing possibility that 
correlations and entanglement may exhibit a completely 
different behavior. Indeed, for gapped spin- 1 sy stems an 
example for this has already been found [Tj. At the 
AKLT [13 point (/3 = -1/3) the entanglement length 
of the ground state diverges, although the correlations 
decrease exponentially. Since the AKLT model is closely 
related to the Heisenberg AF, one might expect that both 
systems show qualitatively the same behavior. 



B. Numerical study 

Let us now have a closer look on the spin-1 Heisenberg 
AF. Since in this case we cannot resort to an analytical 
solution as for the AKLT, we have to rely on numerical 
methods. 

We start our analysis by performing exact diagonaliza- 
tions for an open chain of up to 10 sites. At the endpoints 
we couple to S' = 1/2 spins thus making sure that the sys- 
tem is in the singlet ground state. We are interested in 
the LE between the endpoints of the chain. Since the end 
spins are represented by qubits we can still refer to the 
LE as measured by the concurrence: jv- numer- 
ical analysis shows that the optimal measurement basis 
is given by the same local unitary transformation: 



U ■ 



1 

71 



1 o 1 
o \/2 o 
-1 O 1 



(34) 



as for the AKLT. This strategy produces a maximally 
entangled state between the end spins (if ^ = !)■ This 
surprising result can be understood from the analytical 
study of the AKLT model in [l^. In a singlet valence 
bond picture measurements in the basis (|34|) can be in- 
terpreted as Bell measurements on a (virtual) spin-1/2 
system, which lead to entanglement swapping. 

Let US now investigate, whether this effect depends on 
the choice of the boundary condition or the number of 
sites. For this purpose we apply our numerical method 
outlined in Sect. IV. The LE in its variant Lfj'' (see (g))) 
is computed for a chain with PBC and a large number of 
sites {N = 80). 
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FIG. 7: Calculation of (diamonds) for the ground state 

of the antiferromagnetic spin-1 Heisenberg chain with A'^ = 80 
sites as a function of the spin distance n. For comparison 
we plot the correlation function Q^^'^" computed directly 
from the MPS (squares) and using Monte Carlo (circles). 
Numerical parameters (see Sect. IV): D = 16, MC sweeps 
M = 20, 000. 
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In Fig. [7|we see that the LE saturates at a finite value 
0.960 ± 0.003 for large n, whereas the correla- 
tions decrease exponentially. This demonstrates that the 
ground state of the antiferromagnetic spin-1 Heisenberg 
chain could be used to distribute EPR-like entanglement 
over arbitrary distances by performing focal operations 
on the intermediate spins. As mentioned in Sect. II. D 
this result might be particularly interesting in the con- 
text of quantum repeaters. 

Let US now come back to Haldane's result for the Heisen- 
berg AF stated in the beginning of this section. Our 
numerical study of the spin-1 case might give a first indi- 
cation that, unlike the correlation length, the entangle- 
ment length is infinite for both half-integer and integer 
spins. 



C. Hidden order and string order parameter 

Our numerical results show that the ground state of the 
spin-1 Heisenberg AF exhibits long range order in terms 
of the localizable entanglement. We have also seen that 
this long range order is not reflected in the behavior of 
the two-particle correlation functions. However, one can 
define a multiparticle correlation function, the so called 
string order correlation function 52], which detects this 
hidden order in the ground state. The string order cor- 
relation function has been argued to be of topological 
nature and is defined as 



with Rk = exp [íttS^ 
rameter, ^^o ■= lini„- 



{Sl[(g> 



'i-\-n\ 



(35) 



A non-vanishing string order pa- 
indicates the presence of 



long range (hidden) order. As an obvious generalization 
of the string order correlation function to arbitrary mod- 
els, let US define a connected version in a variational way. 
Consider the set of ali observables {O} with bounded 
spectrum — 1 < O < 1. We define the connected string 
order correlation function Qcgg"' (and the related param- 
eter ^cso) for a given translational invariant state as 



Q 



cso 



max {Ol 

-t<ÔuÔ2<t 



(36) 



Here {A1A2 . . . An)c denotes the connected n-point corre- 
lation function, which can be defined in a recursive way: 

(Al) = (37) 
(Ai^2> = {Ai)M2)c + {AiA2)c, (38) 
= {Ai),{A2)c{A3)c + {A,),{A2As)cm 

+ {A2)c{AiA3)c + {A3)c{AiA2)c 



Note that the connected part assures that Q^so" mea- 
sures a nonlocal correlation, and that the string order 
parameter of the AKLT-ground state is indeed recovered 
by this definition. 

It has been verified numerically that is finite for 
the ground state of the spin-1 Heisenberg AF. This fact 
can rigorously be proven for the related AKLT-ground 
state |23|- For this state it was further shown that the 
LE saturates as well with the spin distance n ^ . Hence 
one might expect a connection between the existence of 
long range order in the entanglement and long range or- 
der indicated by the string order parameter. However, 
one can find examples for which this connection does not 
hold. 

For instance, in Til it has been shown that already an 
infinitesimal deformation of the AKLT model leads to an 
exponentially decreasing LE, whereas Çcso stays finite. 
On the othcr hand, ground states exist that exhibit a di- 
verging entanglement length but vanishing ^cso- A simple 
example can be found in the class of MPS H20|) defined on 
qubits {d = 2) and with qubit bonds {D = 2). Note that 
ali these MPS are guaranteed to be ground states of some 
local Hamiltonians. Furthermore, for MPS with qubit 
bonds the string order parameter and the LE can easily 
be computed analytically (see and Appendix C). In 
particular, let us study the state defined by 



A' 



A' 



il. 



(40) 



- so 



The entanglement length can easily be proven to be infi- 
nite. A necessary condition for ^cso to be nonzero is that 
there exists a unitary operator O2 for which the largest 
eigenvalue of Eq has the same magnitude as the max- 
imal eigenvalue of (Appendix C). For the example 
given, this is impossible, hence providing an example of 
a ground state with a diverging entanglement length but 
no long range hidden order. 



13 



VII. MIXED STATES 

In this section we apply the conccpt of LE to char- 
acterize the entanglement of multipartite mixed states. 
Note that the definition of LE Q already includes the 
possibiUty of having a mixed state p. This imphes that 
the states pj^} — {{s\\p\{s\) after the measurements are 
also mixed. Hence we refer in the foUowing to the LE as 
measured by the negativity: i^. In order to provide a 
tool for the computation of Lf^ let us now generahze the 
numerical method, outhned in Sect. IV, to mixed states. 

A. Numerical method 

The key point is to find a representation of a mixed 
state in terms of low dimensional matrices Ak , analogous 
to the MPS l|2()|l. This problem has been considered re- 
cently in There the concept of MPS is generahzed 
to matrix product density operators (MPDO), which are 
defined as 

d 

p= Tr(Mf '"í . . . ) 

Si,Sj,...,Sjv,S^ = l 

x|si,...,SAr)(s;,...,S^|, (41) 

where M^'"^'' are D\ x Z?^^]^ matrices. They can be 
decomposed as 

M'/ = j2At''®iAi'y- (42) 

a=l 

The state p can be purified into a MPS by including 
ancilla states {!&*;)} of dimension dk.: 

si,...,SN ai,...,aN \k—l / 

_ (43) 
In |4p| a method has been introduced that allows one to 
determine the matrices Ak iteratively for a given Hamil- 
tonian and temper ature. 

Starting from the mixed state in MPDO representation 
(|41|l the LE can be computed along the same hnes, using 
the MC method, as for pure states. However, we note 
that the matrices Mk in 141|l have dimension x D^, 
compared to the DxD matrices Ak for pure states. Thus 
in the case of PBC the computation time for the MC part 
scales at least with d for mixed states, but only with 

for pure states. 
Let US also comment briefly on the optimal measurement 
basis. In contrast to pure states, the optimal basis for 
mixed states can no longer be deduced directly from the 
matrices Mk- Alternatively a good guess for the best 
measurement strategy can be found by exact diagonal- 
ization of the Hamiltonian for small N, foUowed by nu- 
merical optimization of LE. 



B. Example: AKLT model 

Next we apply this numerical method to a specific ex- 
ample. Interesting candidates can be found in gapped 
spin-1 Systems, like the AKLT model or the Heisenberg 
AF (see (|32l)- We have seen that the ground state of 
these two models exhibits an infinite entanglement length 
[lll | . Thus the natural question arises, to what extent this 
feature holds for small but finite temperatures and how 
the entanglement length scales with temperature. 
For our numerical study we choose the AKLT model. 
The reason is simply that, here, mixed states can be ap- 
proximated rather well by matrices with dimensions as 
small as D « 10, even for very small temperatures. This 
makes the computation of LE much more eíRcient com- 
pared to the Heisenberg AF. Efficiency is also the reason 
for choosing OBC. 

Let US now discuss the optimal measurement strategy. 
First of ali we point out, that with OBC the ground 
state of the AKLT is four-fold degenerate 0| . Thus for 
T O the density matrix is an equal mixture of these 
four states, which strongly reduces the LE compared to 
e.g. the singlet ground state studied in 0|. On the 
other hand it is known that the degeneracy results only 
from the cnd spins of the chain. Thus, one can strongly 
reduce this boundary effect by choosing the two spins, 
i and i + n, to be far away from the boundaries. For 
this situation we found that the optimal measurement 
scheme for LE is given by measurements in the í7-basis 
(|34|l on the spins between sites i and i + n, and in the 
standard basis everywhere else. This result is not yery 
surprising in terms of the valence bond picture in [llj. 
Entanglement swapping is only needed between the two 
spins of interest, whereas the effect of the degeneracy can 
be minimized by measuring the outer spins in the stan- 
dard basis. Using this strategy we plot in Fig. IHlthc LE, 
as given by Lf^j,)^^, depending on the spin distance n for 
various temperatures T. The temperatures are chosen to 
be of the order of the enery gap |5^] . 

The data indicates an exponential increase of the en- 
tanglement length ~ e"/"^ with a « 0.8, thus lead- 
ing smoothly to an infinite entanglement length at zero 
temperature. This behavior is not unexpected for a ID 
system from the perspective of the Mermin- Wagner the- 
orem j^. However, it is not clear whether this theorem 
is really applicable to phase transitions in terms of LE. 
It is more inspiring to treat this problem on the basis 
of projected entangled-pair states -òG] . In this picture a 
finite temperature phase transition for LE could possibly 
occur for two or more dimensions. We note that recently 
such a transition has been shown to exist for 3D cluster 
states 14]. 

Finally we would like to point out that, although the en- 
tanglement length of the AKLT model is finite for T > O, 
it can still be considerably large for sufRciently low tem- 
peratures T < 0.2. Thus for practical purposes this sys- 
tem might still be useful, e.g. for quantum repeater se- 
tups. 
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FIG. 8: Calculation of the LE as given by Li ij^n for ^he 
AKLT model as a function of the spin distance n and for var- 
ious temperaturas T. We have chosen a chain with OBC and 
A'^ = 50 sites. We note that almost identical data can be ob- 
tained for N = 2G, indicating that our resuhs are aheady close 
to the thermodynamic hmit. Numerical parameters (see also 
Sect. IV): matrix dimension D = 10, MC sweeps M = 5, 000. 



VIII. CONCLUSION 

In this article we have presented a detailed discussion 
of the idea and basic properties of LE, as defined in 10]. 
Apart from that we also found new results: The central 
finding of Il0l| , namely the lower bound of LE in terms of 
connected correlation functions, has been generalized to 
pure qutrit states. Moreover we have proposed a numer- 
ical method, applicable to both pure and mixed states, 
that allows to calculate the LE efhciently even for large 
ID spin Systems. In future work we are planning to gen- 
eralize this numerical scheme to two and higher dimen- 
sional spin systems based on the method [5a|. We have 
applied our numerical scheme to study the LE of various 
spin models. The results can be summarized as follows: 
The LE, as well as the entanglement fluctuations, ex- 
hibit characteristic features at a quantum phase transi- 
tion. This result is a direct consequence of the numerical 
observation, that, for ground states of spin- 1/2 systems 
with two-spin interactions, the lower bound is typically 
tight, i.e. the LE is completely characterized by the max- 
imal connected correlation function. However, connected 
correlation functions give only a coarse grained picture 
for the LE. As an example, we have shown for the spin- 
1 Heisenberg AF that the entanglement length diverges, 
whereas the correlation length is finite. This numerical 
result might suggest that, in terms of LE and opposed to 
connected correlation functions, no distinction has to be 
made between the scaling behavior of integer and half- 
integer spin Heisenberg AF's. To confirm this idea, fur- 
ther studies of the integer spin case are desired. 
We further note that preliminary results indicate that the 
entanglement features of the Heisenberg AF might hold 
qualitatively for the whole class of gapped spin-1 models 



defined by (|33|l . This finding would allow one to distin- 
guish the Haldane phase from the dimerized phase based 
on the scaling of LE. Let us also mention in this con- 
text that gapped spin-1 systems also appear to be ideal 
candidates to look for quantum phase transitions, being 
detected solely by the entanglement length and not the 
correlation length (like for the generalized AKLT model 
in .11]). 

In order to illustrate that our numerical method works 
also for mixed states we computed the LE of the AKLT 
model for finite temperatures. We found that the en- 
tanglement increases exponentially with the inverse tem- 
perature. This smooth behavior indicates the absence 
of a phase transition in terms of LE for the ID case. 
However, we have reason to believe that such a phase 
transitou might occur for the 2D system. 
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APPENDIX A: CONNECTED CORRELATION 
FUNCTIONS AND LOCAL MEASUREMENTS 
ON PURE QUBIT STATES 

Here we present an extended version of the proof in 
[TÕI ] for the foUowing statement: 

Given a (pure or mixed) state p of N qubits with con- 
nected correlation function between the spins i and 
j and directions a, b, then there always exists a basis in 
which one can locally measure the other spins such that 
this correlation does not decrease, on average. 
Proof: 

Formally we have to show that there exists a measure- 
ment M such that: 



\QÍbÍp'')\ 



s 



(Al) 



To this end let us first consider mixed states of three 
qubits. A mixed 3-qubit density operator can be param- 
eterized by four 4x4 blocks 



Pi o- 



(A2) 



Since local unitary operations can be absorbed in p it 
is sufficient to consider the Ql^ correlations. Thus the 
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original correlations are completely determined by the 
diagonal elements of the reduced density operator pi +/32- 
A von Neumann measurement in the basis 

1+) cos(él/2)|0) +sin(6'/2)e''^|l), (A3) 
|-> := -sin(0/2)e-*^|O) +cos(í?/2)|l>, (A4) 

on the third qubit results in the hormitian unnormalized 
2-qubit operators 



X± ■- (±|p|±) = — - — ±008(6») 



(A5) 



±sin(6l) cos{(t)) 



+ sin((;/>) 



with probabilities p± = Tr(X-i-) conditioned on the out- 
come {+} or {— }. From these equations we see that 
the SU(2) transformation on the third qubit can be ac- 
countcd for by a S0(3) rotation of the z-axis, defined by 
the unit vector: 



X := [cos(é'); sin(é') cos(0); sin(6') sin((/i)] 



(A6) 



As noted above we have to consider only the diagonal 
parts of the measurement outcomes which can be 
represented in terms of the column vectors: 



X± 



1 



R 



1 

±x 



(A7) 



where R is the real 4x4 matrix whose columns consist of 
the diagonal elements of the matrices (pi + P2), {pi ~ P2), 
(ct + ct'!'), «(cr-o-l'). 

Provided with these definitions the inequality (jAlj) can 
be written in the form: 

p+|g,,(X+/p+)|+p_|Q,,(X_/p_)| > \Q,,{X++X^)\. 

Inserting 1 — t'í:{X±/p±), this inequality can be trans- 
formed in a bilinear form in x: 



(1 x^)S 



+ 



1 



(1 ^x^)S 



1 



> |4a! 



Here a is the first element of the matrix 

S := R^{(Jy ® (Jy)R 



(A8) 



(A9) 



and f5,Q are defined as 3 x 1 and 3x3 blocks, respec- 
tively. Without loss of generality we can assume that a is 
positive and thus remove the absolute value sign in lASp . 
Some straightforward álgebra yields then the sufficient 
inequality: 



{A + B)x>{) , 



with 



A := 




B:=Q 



PP 



(AIO) 



(Ali) 



where c is such that p± = (1 ± c^í)/2. We now have 
to show that the matrix A + B has at least one positive 
eigenvalue. From the form of A one immediately sees 
that it is positive semidefinite (a > 0). The matrix B 
rcquires more work. First we note that the matrix Uy^Uy 
in l|A9|l has two negative and two positive eigenvalues. 
Assuming nonsingular R it foUows from Sylvester's law of 
inertia 57] that S also has two positive and two negative 
eigenvalues [õ^, and so has the inverse S^^ . Now B 
is the inverse of the Schur complement of a, and hence 
corresponds to a principal 3x3 block of the matrix S~^: 



s- 



B- 



(A12) 



where the entries * are of no interest here. Let us denote 
the eigenvalues of in algebraic increasing order by 
Al . . . A4 and those of B~^ by /zi . . . From the inter- 
lacing properties of eigenvalues of principal blocks [57j |. 
we obtain the foUowing relation: 

Al < Aii < A2 < Ai2 < A3 < Ai3 < A4. (A13) 

Knowing that A3 > O we deduce that posesses at 
least one positive eigenvalue, and so does B. The exis- 
tence of one positive eigenvalue in A + B ensures that one 
can always find a measurement direction x such that the 
inequality (|A10() is fulfilled. We have proven the theorem 
for a mixed three qubit state. However, this result can 
immediately be extended to arbitrary N . To see this let 
US consider e.g. the correlation \Qzz{X^/p+)\ for one of 
the measurement outcomes on the third qubit. The two 
qubit state Xj^ can be expanded in a basis corrcsponding 
to a measurement of the fourth qubit {X+ = + Y^) . 
The theorem can now be applied with rcspcct to the 
States Y± and so forth, completing the proof. □ 

Note that the proof is constructive and allows to deter- 
mine a measurement strategy that would at least achieve 
the bound reported. 

Let US now show that the above result can also be 
generalized to a setup where the spins i and j can have 
any dimension, but the measurements are still performed 
on qubits. To be more specific we consider the operator 
S\ ® S'g acting on a bipartite state pij of arbitrary di- 
mension. Since local unitary transformations can always 
be absorbed in the definition of pij we can choose S\ and 
5;^ to be diagonal. The correlation function can then be 
written in the bilinear form: 



O'- 



tr[p,j {S'a (S> S'g)] - tr[/5,j {S\ ® l)]tr[p,j (1 ® S^)] 



— -X (a 1 — 1 a ) {g) (6 1 — 1 b )x , 



(A14) 



where the column vectors x, a and b are representing 
the diagonal elements of the matrices pij , S\ and Sg , 
and 1 is a column vector with ali ones. Thus the matrix 
Z := (a l-'" — 1 a^) (g) (ò 1"^ — 1 tí^) replaces the matrix 
CTy (g) CTy in the definition (|A9|) . Z is the tensor product of 
two antisymmetric matrices of rank two and therefore has 
two positive and two negative eigenvalues. This property 
is sufficient to fulfiU the inequality (jA10|) . 
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APPENDIX B: MAXIMUM CORRELATION 
FUNCTIONS AND ENTANGLEMENT OF PURE 
BIPARTITE STATES 

1. Two-qubit states 

For an arbitrary two-qubit state p we want to maximize 
the correlation function: 

Q^AB = tr[p( S\ ® S%)] - tr[p( S\ ® l)]tr[p(l ® S^)]. (BI) 

For qubits we can parametrize Sa and Sb by the three- 
dimensional unit vectors a, b: 



Here we have defined the 3x3 matrices: 



Sa = í? ■ a , 
Sb = a - b , 



(B2) 
(B3) 



where a — {ax Uz)- The correlation can then be writ- 
ten in the form: 



{AB 



a/3 



Qaf3 b 



'/3 



Q b {a,P = X, y, z) 



(B4) 

The matrix elements Qap of the 3x3 matrix Q are defined 
by HB1|I with Sa = Ca, Sb = cfp. Clearly the maximum 
value for Qab is given by the largest singular value of 
the matrix Q. 

For a pure state p = the matrix Q can 

be computed using the Schmidt decomposition j?/;) — 
Tlii ® \^b) with Ai > 0. Note that local unitary 

transformations can always be absorbed in the definition 
of Sa and Sb- In this representation the matrix Q is 
diagonal and one can show that the maximum value is 
given by Qxx = 2A1A2. It can easily be checked that 
this expression is equal to the concurrence C as defined 
in ((SJ- Thus we have shown that the entanglement of a 
pure two-qubit state \'ip) as measured by the concurrence 
C is equal to the maximum correlation function: 



max((3^B(|^))) 



^xim = cm) = 2x1X2 . (B5) 



This relation is central for establishing the strong connec- 
tion between classical and quantum correlations in pure 
multipartite qubit states (see Sect. III). 



2. Two-qutrit states 

As in the qubit case we consider correlations of the 
form ljBl|) . In generalization to the usual spin-1 opera- 
tors we want to maximize with respect to the bounded 
operators —1 < Sa, Sb < 1- In the forthcoming discus- 



sion we will only consider pure states. In Schmidt de- 
composition we have {tp) = J2'i=i N) ^ N) with Ai > O 
and Af -I- A| -I- A§ = 1. Let us begin with rewritting the 
correlation function: 



Qab = í^aÍSa (pi - PD) 
= ítbÍSb ÍP2 - aD) 



(B6) 
(B7) 



Pi 
D 



= trB(l (E>Sb p)^ D^SlD^ , (B8) 
= ívaÍSa ®lp)= DÍSaD^ , (B9) 
= diag{XlXlXl) , (BIO) 



and the scalars a = tr(p2), P = tr(pi). Wc further in- 
troduce the eigenvalue decomposition pi — (3D ~ U EW . 
Note that the diagonal matrix E has zero trace and thus 
has at least one negative entry. Now one immediately 
sees from (IB6|) that Qab is maximized if Sa has the 
same eigenvectors as pi — PD and if its eigenvalues are 
given by the sign of the matrix E. Hence we can formu- 
late the foUowing relations that hold for the maximum 
correlation function: 



(í) 5a = [/ sign(E) Ut 
(ii) [Sa,Pi-PD]=0 
(Ui) [Sb,P2 — ctD] — O . 



with tr(E) = O (Bll) 
(B12) 
(B13) 



These conditions lead to the simple comniutator rela- 
tion [D,M] = O where M := Ui\E\ - a E)W . Since 
AI commutes with the diagonal matrix D it has to be 
diagonal. Trivially this is fulfilled for diagonal U im- 
plying also a diagonal operator Sa- A nondiagonal U 
is only possible if the matrix {\E\ — a E) is degenerate. 
Hence we arrived at the surprising result that the op- 
erator 5*^1 that maximizes the correlation Qab is either 
diagonal Sa = diag(l, — 1, — 1) or can be parameterized 
in the form: 



Sa^ 



1 O O 

o cos(6') sin(6')e-*'^ 
O sin(6')e*'^ - cos(6l) 



(B14) 



As Qab is symmetric in A and B an equivalent ex- 
pression with rotation angles 9' and çi' holds for the 
operator Sb- From this we can deduce the relations 
a = Xl + cos{e'){Xl - Xl) and p = Xl + cos{e){Xl - A^). 

The required degeneracy of the matrix (\E\ — a E) 
puts a constraint on a (or the optimal rotation angle 9') 
as a function of /3: a = F{P) (for simplicity we do not 
specify the function F here). Due to symmetry it also 
holds P = F{a)- Clearly a fixpoint of the maximization 
procedure is given by the symmetric solution a = P (or 
— 9'). However, there also exists an asymmetric solu- 
tion. In order to obtain nice analytical expressions for 
these Solutions it is more conv enient to parameterize the 
function Qab using the form (|B14|) and then maximize 
with respect to the rotation angles O and 6': 

Qab = Aj + cos 6* cos 6''(A2 + A3) + 2 sin sin &'A2A3 

- (A? + (A^-A^)cose)(A? + (AÍ-AÍ)cos0') (B15) 

Here we made the choice (f> — —4>', which maximizes 
Qab- We further note that in this expression the role of 
the Schmidt coefficient Ai is special, which results from 
the ordering Ai < A2 < A3. For the symmetric case 
{9 — 9') we obtain the optimal rotation angle: 



cos 



'opt 



) = 



-^1(^2 ^ ^3) 



(A2-A3)2-(Ai-Ai) 



(B16) 
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which yields the maximum correlation function: 



~.sym _ 
!AB — 



2A2A3 — 



(B17) 



Notice that for Ai = O this reduces to the qubit solution 

Q 



AB 



2A2A3. As for the symmetric case Qab Í|ljl5|) 
is quadratic in cosiO) the maximum can also be reached 
at the boundaries cos{9) = ±1. This leads to diagonal 
operatators 5*^1 = Sb and the maximum correlation is 
given by: 



^AB 



(B18) 



The asymmetric solution can also be worked out, but is 
difhcult to cast in a nice analytical form. For our purpose 
it is enough to establish the foUowing rclation for the 
optimal rotation angles: 

a {cos{e) + cos{e')) = -b{l + cos{e) cos{9') (B19) 
a ■■=XÍ + \l-{\l-Xlf (B20) 



b ■.= xl{xl-xl) 



(B21) 



It can easily be verified that a, 6 > O and a > b. Insert- 
ing l|B19|l in Qab í|B15|) it foUows that the asymmetric 
solution can be upper bounded by: 



Wab S - 



2A9A.':í 



(B22) 



Straightforward analysis shows that the two-qubit limit 
(Al = 0) yields an upper bound for the maximum corre- 
lation in the two-qutrit case: 

Qab ■— inax(y^3 tQab ^Qab i - ^a2A3|ai=o- 

(B23) 

The maximum correlation function decreases if the num- 
ber of non-zero Schmidt coefficients increases. Thus 
Qab^ cannot be used for measuring entanglement as in 
the qubit case. The entropy of entanglement E{\ tp)) [25l| . 
on the contrary, increases with the number of non-zero 
Schmidt coefhcients. A fact that foUows directly from 
the concavity property oí Edip)). Hence the Ai = O case 
E{\ = /(2A2A3|ai=o), with / being the convex func- 
tion lO, yields a lower bound on the entropy of entan- 
glement. From this it foUows that the entanglement of a 
pure two-qutrit state is lower bounded by the maximum 
correlation function: 



Ei\i^))>f{QTBl 



(B24) 



APPENDIX C: ANALYTICAL CALCULATION 
OF THE STRING ORDER PARAMETER AND 
THE LE FOR MATRIX PRODUCT STATES 

We consider a MPS ^ with qubit bonds {D = 2). 
In the case of OBC and qubits at the endpoints 



(i = O, iV -I- 1) this (unnormalized) MPS state can be 
written in the form: 



a,ii...ÍN,0 



.A'-P |a)|^l...^^)|/3), (Cl) 



where a and ò are two dimensional row and column vec- 
tors, respectively, and a,/3 e {0,1}. We are interested 
in the string order parameter (|35|l between the endspins. 
Using expression l(?7jl for calculating expection values of 
MPS, we can write: 



0,JV+1 

so 



í;^_(í;«)^< 



(C2) 



In the limit of large N and diagonalizable E^i [E^ ) only 
the maximum eigenvalue Xn {X^) will survive: 



S.SO := lim Qso 



0,A^+1 



{Êlr^){ltÊ\) \Xt 



N 



(C3) 

where Iq and denote the left and right eigenvectors of 
Eo. 

In the case of the AKLT model and for the basis H34() , we 
have: — iay, = and A'^ — Ox- Hence onc finds 



E^ = (Tx <^crx - 

Er = -CTj, (g) CTj; 



(C4) 
(C5) 



and Xr — X^ ^ 3. Realizing that a" and ír are repre- 
senting unit vectors in the standard basis, we obtain the 
result: Çso — 1- 

Let US now show how to calculate the LE between 
the end points of the chain for states of the form I|C1|I . 
Since the end spins are represented by qubits we can 
use the concurrence (O as entanglement measure, which 
simpliiies the calculation considerably. For the basis 
A4 — the average entanglement can be written 

as [13: 



.M,c _ E.,.....2|det(^'^...A'")| 
ÊliE^fÊl 

Since the determinant factorizes, we obtain 



E, 



M,C 
0,^ + 1 



{Elr^){l^E\) V Al 



EJdet(A*)| 



N 



(C6) 



(C7) 



in the limit of large N . The basis which maximizes L^ '~^ 
is clearly the same basis, which maximizes the expression 
|det(^')|. This problem is equivalent to calculating 
the EoA of the D"^ x state A+A: 

Ea{Á) :=supy |det(A*)| = ir\A^ {ay ® (jy)A\. (C8) 
^ i 

The elements of the (25 + 1) x matrix A are given by 
Hence we found a necessary and sufR- 



,("/5) 



cient condition for long range order in the entanglement 
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(i.e. non-vanishing i[J^jv+i ^'^^ ^ ~^ °°)'- expres- 
sion EaÍA) has to be equal to the largest eigenvalue, X-^, 
of the matrix For the AKLT model, one can eas- 
ily check that this condition is indeed fulfiUed, and that 



The ground state of the AKLT thus exhibits long range 
order both in terms of the LE and the string order pa- 
rameter. 
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